/*******************************************************************************
*   			Impacts of a Large-Scale Parenting Program: 				   *
*					Experimental Evidence from Chile						   *
********************************************************************************


	REQUIRES:	"${dt_final}panel_RR_indexes.dta"
	CREATES:	${out_tables}/T7A_ITT_child_individual.xls
				
	WRITEN BY:  Italo Lopez  [italolop@usc.edu]

********************************************************************************
	Prepare data
*******************************************************************************/

	* Load data
	use "${dt_final}panel_RR_indexes.dta", clear


	* Establish macros for regressions

/*******************************************************************************
*   Set macros for regressions
*******************************************************************************/
	
	* Outcomes
	*-----------------------------------------------------------
	local childoutcomes  "zDCCSPunt zptevir_irtscore2 zt_score zcbcl_t_e zcbcl_t_i zibatt_sspi zibatt_ssai zibatt_sssr zas_raw_score" 
	local behaviors "home_score se_score disc_score" 
	local beliefs "zpstyle1 zpstyle2 zpstyle3 zPSCS zPSSS_family zPSSS_friends zPSSS_others zPSI_Distress zCESD zpacotis"	
	local child "zptevir_irtscore2 zcog_index zsoc_index" 

	* Controls
	*-----------------------------------------------------------

	macro  def Xvar1 "i.age_year i.gender"  	
	macro  def Xvar2 "i.age_year i.gender inc_qtaut_old i.pc_edu_mdsfin hh_mem hh_tipo"  	
	macro  def Xvar2old "i.gender inc_qtaut_old i.pc_edu_mdsfin_old hh_mem hh_tipo"  

	
	* Parental characteristics: WAIS and BFI and imputed wais 
	*-----------------------------------------------------------

	macro def parent_cha4 "BFI_ext BFI_ope BFI_con BFI_neu BFI_agr wais2 flag_wais"
	

/*******************************************************************************
	Table A7: ITT estimates of child individual measures
*******************************************************************************/



	local measures "zptevir_irtscore2 zDCCSPunt zas_raw_score zibatt_sspi zibatt_ssai zibatt_sssr zcbcl_t_t zcbcl_t_e zcbcl_t_i zcbcl_p_1 zcbcl_p_2 zcbcl_p_3 zcbcl_p_4 zcbcl_p_5 zcbcl_p_6 zcbcl_p_7"
	local list1 "measures" 
	#delimit;
	foreach f of local list1 {;

			foreach var of local `f' {;
			
				di "`var'";
					
				* PANEL 1: ITT with Age and Gender ;
					xi: xtreg `var' i.TIPO $Xvar1,  fe i(CENTRO_SALUD);
					estimates store `var'_itt1;
					test _ITIPO_2==_ITIPO_3;
					estadd scalar pb=r(p);				

				* PANEL 2: ITT with Age and Gender and parental characteristics ;
					xi: xtreg `var' i.TIPO $Xvar2 $parent_cha4,  fe i(CENTRO_SALUD);
					estimates store `var'_itt2;
					test _ITIPO_2==_ITIPO_3;
					estadd scalar pb=r(p);		
					
				* PANEL 3: ITT plus baseline outcome ;
					xi: xtreg `var' i.TIPO $Xvar2 $parent_cha4 zchild_index_old,  fe i(CENTRO_SALUD);
					estimates store `var'_itt3;
					test _ITIPO_2==_ITIPO_3;
					estadd scalar pb=r(p);				
					
								};
			
	};


	#delimit cr
	***  Tables ITT Outcomes

	
		esttab zptevir_irtscore2_itt1 zDCCSPunt_itt1 zas_raw_score_itt1       /*
		*/  zcbcl_t_t_itt1  /*
		*/ zcbcl_t_e_itt1 zcbcl_t_i_itt1 zcbcl_p_1_itt1 zcbcl_p_2_itt1        /*
		*/ zcbcl_p_3_itt1 zcbcl_p_4_itt1 zcbcl_p_5_itt1 zcbcl_p_6_itt1        /*
		*/ zcbcl_p_7_itt1   zibatt_sssr_itt1 zibatt_ssai_itt1 zibatt_sspi_itt1    	  		   								  /*
		*/ using "${out_tables}/appendix/T7A_ITT_child_individual.csv", replace f		  /* 
		*/ keep(_ITIPO_2 _ITIPO_3) 											  /*
		*/ cells("b(fmt(3)star)" "se(fmt(3)par)") 							  /*
		*/ star(* 0.10 ** 0.05 *** 0.01) 			legend					  /*
		*/ stats(r2 pb N, fmt(2 3 0) labels("R-squared" "P value for B=I"  "Observations")) 	 /*
		*/ mtitles("TEVI-R (language)""DCCS (cognitive flexibility)" 		  /*
		*/ "Leiter-R (sustained attention)" "CBCL total (behavior)" 		  /*
		*/ "CBCL externalization" "CBCL internalization" "CBCL Emotional"	  /*
		*/ "CBCL Anxious""CBCL Somatic" "CBCL Withdrawn" "CBCL Sleep problems" /*
		*/ "CBCL Attention problems" "CBCL Aggressive" "Battelle Social Role" /*
		*/ "Battelle Interaction w/Adults" "Battelle Interaction w/Peers" )   /*
		*/ title("Child Individual Outcomes (1): (age/gender)") 	varlabels(_ITIPO_2 "NEP-B" _ITIPO_3 "NEP-I")
		

		esttab zptevir_irtscore2_itt2 zDCCSPunt_itt2 zas_raw_score_itt2		  /*
		*/  zcbcl_t_t_itt2  /*
		*/ zcbcl_t_e_itt2 zcbcl_t_i_itt2 zcbcl_p_1_itt2 zcbcl_p_2_itt2 		  /*
		*/ zcbcl_p_3_itt2 zcbcl_p_4_itt2 zcbcl_p_5_itt2 zcbcl_p_6_itt2		  /*
		*/ zcbcl_p_7_itt2 	zibatt_sssr_itt2 zibatt_ssai_itt2 zibatt_sspi_itt2												  /*
		*/ using "${out_tables}/appendix/T7A_ITT_child_individual.csv", append f		  /* 
		*/ keep(_ITIPO_2 _ITIPO_3) 											  /*
		*/ cells("b(fmt(3)star)" "se(fmt(3)par)") 							  /*
		*/ star(* 0.10 ** 0.05 *** 0.01) 			legend					  /*
		*/	stats(r2 pb N, fmt(2 3 0) labels("R-squared" "P value for B=I"  "Observations")) 		 							 /*
		*/ mtitles("TEVI-R (language)""DCCS (cognitive flexibility)" 		  /*
		*/ "Leiter-R (sustained attention)" "CBCL total (behavior)" 		  /*
		*/ "CBCL externalization" "CBCL internalization" "CBCL Emotional"	  /*
		*/ "CBCL Anxious""CBCL Somatic" "CBCL Withdrawn" "CBCL Sleep problems" /*
		*/ "CBCL Attention problems" "CBCL Aggressive" "Battelle Social Role" /*
		*/ "Battelle Interaction w/Adults" "Battelle Interaction w/Peers" )   /*
		*/ title("Child Individual Outcomes (2): + caregiver's education, IQ and personality traits; hh'ld income, hh'ld size") 	varlabels(_ITIPO_2 "NEP-B" _ITIPO_3 "NEP-I")

		esttab zptevir_irtscore2_itt3 zDCCSPunt_itt3 zas_raw_score_itt3 	  /*
		*/  zcbcl_t_t_itt3  /*
		*/ zcbcl_t_e_itt3 zcbcl_t_i_itt3 zcbcl_p_1_itt3 zcbcl_p_2_itt3 		  /*
		*/ zcbcl_p_3_itt3 zcbcl_p_4_itt3 zcbcl_p_5_itt3 zcbcl_p_6_itt3 		  /*
		*/ zcbcl_p_7_itt3 	zibatt_sssr_itt3 zibatt_ssai_itt3 zibatt_sspi_itt3												  /*
		*/ using "${out_tables}/appendix/T7A_ITT_child_individual.csv", append f		  /* 
		*/ keep(_ITIPO_2 _ITIPO_3) 											  /*
		*/ cells("b(fmt(3)star)" "se(fmt(3)par)") 							  /*
		*/ star(* 0.10 ** 0.05 *** 0.01) 			legend						  /*
		*/ stats(r2 pb N , fmt(2 3 0)  /*
		*/ labels("R-squared" "P value for B=I" "Observations" )) 	  /*
		*/ mtitles("TEVI-R (language)""DCCS (cognitive flexibility)" 		  /*
		*/ "Leiter-R (sustained attention)" "CBCL total (behavior)" 		  /*
		*/ "CBCL externalization" "CBCL internalization" "CBCL Emotional"	  /*
		*/ "CBCL Anxious""CBCL Somatic" "CBCL Withdrawn" "CBCL Sleep problems" /*
		*/ "CBCL Attention problems" "CBCL Aggressive" "Battelle Social Role" /*
		*/ "Battelle Interaction w/Adults" "Battelle Interaction w/Peers" )   /*
		*/ title("Child Individual Outcomes (3): + baseline outcomes") 	varlabels(_ITIPO_2 "NEP-B" _ITIPO_3 "NEP-I")
		
	
